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A nitrogen-vacancy (NV) center in diamond is a promising sensor for nanoscale mag¬ 
netic sensing. Here we report electron spin resonance (ESR) spectroscopy using a 
single NV center in diamond. First, using a 230 GHz ESR spectrometer, we per¬ 
formed ensemble ESR of a type-lb sample crystal and identified a substitutional sin¬ 
gle nitrogen impurity as a major paramagnetic center in the sample crystal. Then, 
we carried out free-induction decay and spin echo measurements of the single NV 
center to study static and dynamic properties of nanoscale bath spins surrounding 
the NV center. We also measured ESR spectrum of the bath spins using double 
electron-electron resonance spectroscopy with the single NV center. The spectrum 
analysis of the NV-based ESR measurement identified that the detected spins are 
the nitrogen impurity spins. The experiment was also performed with several other 
single NV centers in the diamond sample and demonstrated that the properties of 
the bath spins are unique to the NV centers indicating the probe of spins in the 
microscopic volume using NV-based ESR. Finally, we discussed the number of spins 
detected by the NV-based ESR spectroscopy. By comparing the experimental result 
with simulation, we estimated the number of the detected spins to be < 50 spins. 
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I. INTRODUCTION 


A nitrogen-vacancy (NV) center is a paramagnetic defect center in diamond. A NV center 
is a great testbed to investigate quantum physics because of its unique electronic, spin, and 
optical properties including its stable fluorescence (FL) signals,-i long decoherence tirne,- - - 
and capability to initialize the spin states of NV centers by applying optical excitation and 
to readout the states by measuring the FL intensity.- A NV center is also a promising 
magnetic sensor- - — because its extreme sensitivity to the surrounding electronical and 
nuclear spins.-d Spin sensitivity of electron spin resonance (ESR) spectroscopy is drasti¬ 
cally improved using NV centers. ESR spectrum of small ensemble electron spins has been 
measured using a double electron-electron resonance (DEER) spectroscopy with single NV 
centers. 21- — ESR detection of a single electron spin has also been demonstrated using the 
DEER technique.— - — In addition, bio-compatibility and excellent chemical/mechanical sta¬ 
bility of diamond makes a NV center suitable for applications of nanoscale magnetic sensing 
and magnetic resonance spectroscopy in biological systems.— - — 

In this article, we discuss ESR spectroscopy of a small ensemble of electron spins using 
single NV centers in diamond. Although the state-of-the-art of NV-based ESR spectroscopy 
is the detection of a single electron spin, the small ensemble measurement is often advan¬ 
tageous for applications of NV-based ESR because of less sophisticated sample preparation 
(e.g. higher density of the target spins which increases the coupling to the NV) and more 
sensitive detection (e.g. more pronounced changes in the coherence (T 2 ) and population (T\) 
decays), and still enables to probe nanoscale local environments which may be different from 
the bulk properties. On the other hand, it is challenging to estimate the detected number 
of spins from the small ensemble measurement. 

In the investigation, our experiment is performed with a type-lb diamond crystal at 
room temperature. Using a 230 GHz ESR spectrometer, we first perform ensemble ESR of 
the sample crystal and End that a major paramagnetic impurity in the sample crystal is 
a substitutional single nitrogen center (N spin; also known as PI center). Next, we detect 
a single NV center using FL autocorrelation and optically detected magnetic resonance 
(ODMR) measurements. We then carry out Rabi, free-induction decay (FID) and spin echo 
(SE) measurements of the single NV center to study static and dynamic properties of bath 
spins surrounding the NV center. We also employ DEER spectroscopy to measure ESR 


2 



(a) 



1 mm 



FIG. 1. (a) Picture of the type-lb diamond crystal used in the investigation, (b) cw ESR spectrum 

of N spins measured using the 230 GHz ESR spectrometer. The spectrum was obtained by a single 
scan at 0.2 rnT/s with field modulation of 0.03 mT at 20 kHz. 

signals of the bath spins using the NV center. The detected bath spins are identified as 
N spins from the analysis of the observed ESR spectrum. Based on the intensity of the 
observed NV-based ESR signal, the detected magnetic held by the DEER measurement is ~ 
6 /iT (equivalent to the magnetic held caused by a single S — 1/2 spin with the distance of 
~ 7 nm). We also study several single NV centers in the same crystal and find heterogeneity 
in their spin properties which indicates that our ensemble measurement still proves bath 
spins in the microscopic scale. Finally, we introduce a computational simulation method 
for the estimate of the number of the detected spins. By comparing the experimental result 
with simulation, we estimate the number of the detected spins to be < 50 spins. 

II. RESULTS AND DISCUSSION 

We studied a single crystal of high-temperature high-pressure type-lb diamond, which is 
commercially available from Sumitomo electric industries. The size of the diamond crystal 
is 1.5 x 1.5 x 1 mm 3 (see Fig. [T|(a) inset). The concentration of N spins is 10 to 100 ppm, 
corresponding to 4 x 10 15 to 4 x 10 16 N spins existing in diamond. First, using a 230 GHz 
ESR spectrometer,— we measured ensemble ESR of the sample diamond crystal at room 
temperature to characterize its bulk properties where the magnetic field was applied along 
the (lll)-direction of the single crystal diamond in the measurement. As shown in Fig. [0(b), 
continuous-wave (cw) ESR spectroscopy revealed five-line ESR signals corresponding to the 
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N spins. The ESR intensity of the N spins was drastically stronger than the remaining 
signals, which indicates that the number of N spins dominates the spin population in the 
sample. Moreover, no ESR signals from NV centers and other paramagnetic impurities were 
observed in the ESR measurement because of their low concentration in the sample crystal. 
The spin Hamiltonian of N spin is given by 

H n = fi B g N S N ■ B 0 + S N • A N I N - ^g%I N • B 0 + Pf(/f ) 2 , (1) 

where /ig is the Bohr magneton, B 0 is the external magnetic field, and are the re¬ 
values of the N electron spin, and S N and I N are the electronic and nuclear spin operators, 
respectively. A N is the anisotropic hyperhne coupling tensor to 14 N nuclear spin. The 
gyromagnetic ratio of 14 N nuclear spin (fi^g^/h) is 3.077 MHz where h is the Planck 
constant. The last term of the Hamiltonian is the nuclear quadrupole couplings. The 
previous studies measured g^ y = 2.0024, g* = 2.0025, A£ = Ay = 82 MHz, A^ = 114 
MHz and = —5.6 MHz.- 1 ^ 1 ^ As shown in Fig. QJb), the experimental data agree well 
with simulated ESR spectrum using Eq. (PQ) and the previously determined parameters. 

Figure [2] shows ODMR measurements of a single NV center in the diamond crystal. 
The ODMR experiment was performed using a home-built confocal microscope system.— 
For microwave excitations, two microwave synthesizers, a power combiner, and a 10 watt 
amplifier were connected to a 20 pm-thin gold wire placed on a surface of the diamond 
sample. First a FL image of the diamond was recorded in order to map out FL signals 
in the diamond crystal (see Fig. [21(a)). After choosing an isolated FL spot, we carried out 
the autocorrelation and cw ODMR measurements in order to identify FL signals from a 
single NV center. As shown in Fig. [21(b), the autocorrelation measurement of the chosen FL 
spot revealed the dip at zero delay which confirmed the FL signals originated from a single 
quantum emitter. In addition, cw ODMR measurement of the selected single NV center 
was performed with application of the external magnetic field (Bo) of 35.7 mT along the 
(111) axis. As shown in Fig. [21(c), the reduction of the FL intensity was observed at the 
microwave frequency of 1.868 GHz corresponding to the ODMR signal of the ms = 0 -H- — 1 
transition of the NV center. Thus, the observations of the autocorrelation and cw ODMR 
signals confirmed the successful identification of the single NV center (denoted as NV 1). 

Next, we performed pulsed ODMR measurements of NV 1. In the pulsed measurements, 
a NV center was first prepared in the ms = 0 state by applying an initialization laser pulse, 
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FIG. 2. ODMR experiment of NV 1. (a) FL image of the diamond crystal. The scanned area 
is 3x3 /irn 2 . The color scheme for FL intensity is shown in the legend. NV 1 is circled, (b) 
the autocorrelation data of NV 1. (c) cw ODMR experiment. The solid line indicates a fit to 
the Lorentzian function, (d) Rabi oscillation experiment at Bq = 35.7 rnT and the microwave 
frequency of 1.868 GHz. Square dots connected with lines indicate the measurement and solid line 
indicates the fit. The oscillation frequency ( fn a bi ) of 7.4 ± 0.2 MHz and the decay time (T^) of 0.77 
± 0.04 n s were obtained from a fit to 1/2[cos(27t fnautp) exp(— (tp/Td) 2 ) + 1]. The inset shows the 
pulse sequence with excitation laser (Exc.), microwave (MW) and FL measurement (FL). (e) FID 
measurement at B$ = 35.7 mT and 1.868 GHz. The length of 7r/2-pulse was 34. FL intensity decay 
was monitored as a function of t. The inset shows the pulse sequence used in the measurement, 
(f) SE experiment at Bq = 35.7 mT and 1.868 GHz. The lengths of tt/2- and 7r-pulse used were 
34 and 68 ns, respectively. FL intensity decay was monitored as a function of 2r. The inset shows 
the pulse sequence used in the measurement, b and tc were extracted by fitting FID and SE data 
with Eqs. m and ((3|). For pulsed ODMR measurements, laser initialization pulse (Init.) of 2 /xs, 
and laser read-out pulse (RO) and FL measurement pulse (Sig.) of 300 ns were used. Also the FL 
intensity was normalized and re-scaled into the ms = 0 state population of NV 1. 

and the microwave pulse sequence was applied for the desired manipulation of the spin state 
of the NV center, then the final spin state was determined by applying a read-out laser pulse 
and measuring the FL intensity. In addition, the FL signal intensity was mapped into the 
population of the NV’s ms = 0 state ( P(ms = 0)) using two references (the maximum and 
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minimum FL intensities corresponding to the ms = 0 and ms = — 1 states, respectively).— 
Moreover, the pulse sequence was averaged on the order of 10 6 -10‘ times to obtain a single 
data point. First, the Rabi oscillation measurement was performed at B 0 = 35.7 mT and 
1.868 GHz, which corresponds to the ms = 0 -H- —1 transition of NV 1. As the microwave 
pulse length (tp) was varied in the measurement, pronounced oscillations of P(ms = 0) was 
observed as shown in Fig. [2])c). By fitting the observed Rabi oscillations to the sinusoidal 
function with the Gaussian decay envelope,— the lengths of n/2- and 7r-pulses for NV 1 were 
determined to be 34 and 68 ns, respectively. Second, FID was measured using the Ramsey 
fringes. The pulse sequence are shown in the inset of Fig. [2(e). As seen in Fig. [2))e), the FL 
intensity of NV 1 was recorded as a function of free evolution time (t) and FID was observed 
in the range of t < 100 ns. Third, we carried out the SE measurement. Figure |2](f) shows the 
FL intensity as a function of free evolution time (2r). As shown in the inset of Fig. [2jf ), the 
applied microwave pulse sequence consists of the conventional spin echo sequence, widely 
used in ESR spectroscopy,— and an additional 7r/2-pulse at the end of the sequence to 
convert the resultant coherence of a NV center into the population of the ms = 0 state.- 
The FID and SE decay in electron spin baths have been successfully described by treating the 
bath as a classical noise held ( B n (t )) where B n (t ) was modeled by the Ornstein-Uhlenbeck 
(O-U) process with the correlation function C(t) = (B n (0)B n (t)) = b 2 exp(— \t\/rc), where 
the spin-bath coupling constant ( b ) and the rate of the spin hip-hop process between the 
bath spins (1 The FID and SE decay due to the O-U process are given by, 

FID(t) = 1 - hi + 2cos(27ivpn)] exp[—(ferc) 2 (— + - l)], (2) 

2 6 r c 

and 

SE(t) = - + - exp [-{br c ) 2 {— - 3 - e - (2r)/rc + 4 e - (2r)/(2rc) )], (3) 

2 2 t c 

where A^ = 2.3 MHz is the hyperhne coupling of NV center.— In the quasi-static limit 
( br c 1) indicating slow bath dynamics, SE(t ) ~ exp[— b 2 t 3 /(12T C )\ = exp[—(f/T 2 ) 3 ] 
where T 2 is the spin decoherence time.——— By fitting both FID and SE data (Fig. [2](e) 
and (f)) simultaneously with Eqs.(|2]) and ([3]), we obtained b and Tc to be 30 ± 4 (rad//is) 
and 144 ± 39 /rs, respectively. The result indicates that the surrounding spin bath is in the 
quasistatic limit (brc = 4320 1), therefore, T 2 (= (12r c /5 2 ) 1//3 ) of NV 1 is 1.2 /rs. Using 

the previous study,— the local concentration of the bath spins around NV 1 is estimated to 
be 20 ppm. 
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FIG. 3. NV-based ESR experiment using DEER spectroscopy, (a) The pulse sequence of the 
DEER spectroscopy. In addition to the spin echo sequence of the NV center (denoted as MW), 
another microwave pulse (denoted as MW2) was applied 12 ns after the 7r-pulse in MW to rotate 
the surrounding spins, (b) The obtained ESR spectrum. In the measurement, r = 500 ns and 
tMW 2 = 90 ns were used. The simulated ESR spectrum is shown in the dotted line. 

Furthermore, we performed NV-based ESR using DEER spectroscopy at Bo = 35.7 mT 
(Imwi =1-868 GHz corresponding to the ms = 0 fi —1 transition). Figure |3](a) shows 
the pulse sequence used in the DEER measurements,—which was adopted from the field 
of ESR.— The DEER sequence consists of the spin echo sequence for a NV center and an 
additional 7r-pulse at different microwave frequency (denoted as MW2 in Fig. [3](a)). When 
the additional 7r-pulse is resonant with surrounding electron spins near the NV center, the 
magnetic moments of the surrounding spins are flipped, and this alters the magnetic dipole 
field experienced by the NV center from the surrounding spins during the second half of the 
spin echo sequence. The alteration causes a shift in the Larmor frequency of the NV center, 
which leads to different phase accumulation of the Larmor precession during the second half 
of the sequence from the first half. As a result, the NV center suffers phase-shift in the 
echo signal and reduction of FL intensity from the original spin echo signal is observed. As 
shown in Fig. [3](b), we monitored the spin echo intensity of NV 1 at a fixed r as a function of 
Imw 2 , and observed clear intensity reductions at five frequencies (/ mw 2 — 0.90, 0.92, 1.01, 
1.10, and 1.12 GHz). In the measurement, r = 500 ns and tMW 2 — 90 ns were chosen to 


7 








(a) (b) (c) 







MW2 Frequency (f MW2 ) (GHz) 


FIG. 4. SE, FID and NV-based ESR data of other three single NV centers: (a)NV 2; (b) NV 3; 
(c) NV 4. 

maximize the NV-based ESR signals. As shown in Fig. [3](b), the resonant frequencies of the 
observed NV-based ESR signals were in a good agreement with the ESR of N spins calculated 
from Eq. (Q]), which confirms the observation of N spin ESR spectrum. The DEER intensity 
(Ideer) represents the change of the SE intensity given by the change in the effective dipolar 
field (Bdi Pj eff)- Assuming that 7r-flip of the N spins is instantaneous (he. the effect of a 
hnite width of the pulse was not considered), Ideer = cos(gNv^B(‘I‘'r)Bdip, e ff/K) where 
gisrv = 2.0028 is the g-value of the NV center,— h is the reduced Planck constant, and 
r = 500 ns in the present case. In addition, including the effects of T 2 decay, the NV-based 
ESR signal in P(ms = 0) is given by, 

Invesr = 1/2(1 + I DEER SE(2t)). (4) 

Using Eq. (|4|), we obtained that the intensity of the NV-based ESR at Jmw 2 = 0.90 GHz 
(Invesr ~ 0.625) corresponds to B d i Pt eff ~ 6 /iT. This magnetic field is equivalent to the 
dipole magnetic field from a S = 1/2 single spin with the NV-spin distance ( d ) of ~7 nm 
(Bdip = ho/(47r)(3n(n • m) — m)/d 3 ) with n//m, |m| = ggsms, 9 = 2 and ms = 1/2). 

In addition to the investigation of NV 1, we have studies other NVs in the same diamond 
crystal. Figure |4] shows the set of FID, SE and DEER data from other three NV centers 
(labeled as NV 2-4). The results from NV 1-4 were different because of their heterogeneous 
nanoscale local environments. As shown in Fig. |4](a)-(c), we found that the obtained tq 
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values were varied widely from 144 /is to 3908 {i s. On the other hand, the variation of b 
= 30-38 rad//is is much smaller than re■ As the result, T 2 of NV 1-4 ranges from 1.2 to 
3.4 {is. We also noticed that all NV 1-4 are in the quasistatic limit (brc 1). As shown 
in Fig. [3](b) and Fig. [J the NV-based ESR spectra of NV 1-4 were also quite different. 
The NV-based ESR spectra from NV 1, 2 and 4 displays five-line ESR signals (Fig. [3(b) 
and Fig. [4](a)(c)) whereas NV 3 only shows three visible signals at 0.381 0.474 and 0.553 
GHz (Fig. 13(b)). Possible reasons of the difference are heterogeneity of the number and the 
spatial configuration of N spins aroimd the NV centers, and an uneven distribution of the 
N spin orientations due to a small number of N bath spins. 

Finally, we analyze the intensity of the observed NV-based ESR signal to estimate the 
number of the detected N spins. Our analysis here focuses on the signal at /mw 2 = 0.90 GHz 
corresponds to the \ms = —1/2, mi = —1) \ms = 1/2, m/ = —1) transition of N spins 
oriented along the (111) direction (Fig. [3(b)). In the analysis, we simulate DEER intensities 
by calculating the effective magnetic dipole held at the NV center from the surrounding N 
spins ( Bdip } eff )• First step of the simulation was to generate a model configuration of the 
NV center and N spins in a diamond lattice, and this was done by placing the NV center 
at the origin of the diamond lattice and assigning the positions of N spins randomly in 
the lattice sites (see Fig. [3(a))- Based on the T 2 measurement (T 2 = 1.2 {is ~ 20 ppm), 
~18,000 N spins were placed in the simulated diamond lattice with a diameter of ~200 nm . 
Next, we randomly chose the orientation of N spins where only a quarter of N spins were 
assigned along the (111) direction and the rest of N spins were assigned along the other three 
directions, i.e., [Ill], [111], and [111]. We then assigned the nuclear spin value (mi) of either 
1, 0, or —1 to all N spins with equal probability. For the simulation, we only considered the 
contribution coming from the N spins oriented along the (111) direction with mi = —1 (The 
signal at Jmw 2 = 0.90 GHz in Fig. [3(b)). Thus, B dipe ff was computed from only 1 /12-th 
of all N spins in the simulated lattice on average (~1500 N spins). Finally we assigned the 
electron spin value (ms) of either 1/2 or —1/2 with equal probability to the N spins. 

Bdip,eff is given by the sum of individual dipolar held from each N spin as, 


n N 


B. 


dip,eff 




dip,i 


(5) 


i= 1 

where B diPti = ^ ./'~ 1> g N {i B m l s is the dipolar held strength of z'-th N spin at the NV 
center, un is the number of the N spins oriented along the (111) direction with m d = —1, po 
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FIG. 5. Simulation of the NV-based ESR signals, (a) Overview of the simulation model. Diameter 
of the simulated lattice is 600a where a = 0.375 nm is the lattice constant of diamond. Red 
sphere denotes the NV center placed at the origin of the lattice. Blue spheres represent randomly 
distributed N spins with 20 ppm concentration. Green spheres are a subset of N spins which are 
in the (111) orientation with mj = —1. N spins at or next to the origin, next to other N spins, 
and overlap to other N spins, were suppressed in the simulation, (b-d) Intensity plots of (n/v, 90 %) 
versus Invesr for T\ n = 1 s, 10 ms, and 100 s, respectively. 10 4 spatial configurations were 
simulated for each T\ n . The color scheme for intensity is shown in the legend. Vertical solid lines 
indicate the intensity of the observed ESR signal, namely Invesr = 0.60-0.65. The insets show the 
occurrence of spatial configurations of the simulated lattice that resulted in 0.60 < Invesr < 0-65 
as a function of (rpv, 90 %)- The number of such configurations were ~1200, ~1300, and ~ 400 for 
Ti n = 1 s, 10 ms, and 100 s, respectively, which translate to 4—13 % out of total 10 4 configurations. 


is the vacuum permeability, r* and are the magnitude and polar angle of the vector that 
connects the NV center and z-th N spin, respectively, and m l s is the electron spin value of 
i-th N spin. The mutual flip-flops within N electron spins was also not considered because of 
a low rate of the spin flip-flop process during the measurement, i.e. re t (see Fig. |2](e)). 
We then calculated the number of N spins (un, 90 %) out of Un to obtain more than 90% of the 
Bdip,eff by successively adding individual dipolar field term from each N spin in descending 
order of magnitude (i.e., when B dipe ff — Y^=i° % B c u P ,] |/1 B dip , eff | < 0.1). In addition, we 
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took into account for the electron and nuclear spin relaxations of N spins during the DEER 
measurement time (~ 100 s with ~ 10 7 averaging) by statistically re-assigning the electron 
and nuclear spin values (roj and mj) according to the electron and nuclear longitudinal 
relaxation times ( T le and Tin), respectively. Although T le has been investigated previously 
(T le ~10 ms),- 1 ^ Ti n has not been reported before to the best of our knowledge. In the 
simulation, we used T\ e = 10 ms and considered Ti n to be 10 ms, 1 s and 100 s, therefore 
the assignment of m/ happens 10 2 times for T\ n = 1 s (10 4 and 1 times for T\ n = 10 ms and 
100 s, respectively) and a total of 10 4 iterations (100 s/Ti e ) occurs during the simulation. 
Then, what we finally computed were the average of 10 4 values of un,9 0 % and I deer (he., 
(n N ,90%) = Eli i n N,90%o and i 1 deer ) = El=i 7 dbbr)i therefore, by rewriting Eq. flD, the 
simulated NV-based ESR signal is given by, 

Invesr = 1/2(1 + (Ideer) S E(2r)). (6) 

In addition, in order to see the spatial configuration dependence on (un,90%) and Invesr, 
we repeated the procedure described above for 10 4 spatial configurations and calculated 10 4 
values of (n Nj90 %) and Invesr- 

Figure [5](b)-(d) shows the simulated result of (un, 90 %) as a function of Invesr (Eq. (j6j)) 
for T\ n = 1 s, 10 ms, and 100 s, respectively, from 10 4 spatial configurations. As shown in 
Fig- EXb) , the number of the N spins that contribute to the NV-based ESR signal (( 77 .^ 90 %)) 
depends on the detected Invesr- We first noted that, for all 10 4 configurations we simulated, 
only a small portion (<60) out of ~1500 N spins contributes to more than 90 % of Bdi P , e ff 
on average. In the case of homogeneously distributed N spins, the 60 spins are located 
in a sphere with ~70 nm diameter. Moreover, when the NV-based ESR intensity is large 
(Invesr ~ 0.5), (un, 90 %) is smaller because Bdi P , e ff in such spatial configurations is only 
dominated by a smaller ensemble of the N spin located in the vicinity of the NV center. In 
addition, there was little observation of ( Ideer) < 0 (i.e., Invesr < 0.5) because a small 
probability exists for the NV center to be coupled with a single or a few N spins with the 
20 ppm N concentration and the resultant (Ideer) is a weighted sum of many oscillatory 
functions with different frequencies. On the other hand, when the NV-based ESR intensity 
is small ((Ideer) ~ 1) (he. Invesr ~ 0.73), (un,90%) is larger because the N spins in 
such spatial configurations spread uniformly and the N spins located farther away from the 
NV also contribute to B d i Pt eff- The experimentally observed NV-based ESR intensity was 
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0.625 ± 0.025 (see Fig. [3](b) ). The inset of Fig. 0T) shows a histogram of the occurrence 
of (tt-tv, 90 %) from the simulation that yielded Invesr = 0.6-0.65. The occurrence was in the 
range of 7 — 42 for T ln = 1 s. Moreover, as shown in Fig. 0(c), (riN$o%) with T ln = 10 ms 
was similar to the result with T ln = 1 s (The occurrence was in the range of 7 — 45). On 
the other hand, in the case of Tin = 100 s, the distribution was slightly different and the 
occurrence was in 2 — 28 as shown in Fig. 0(d). Thus, from the simulation, the number of 
N spins in the present NV-based ESR measurement is estimated to be < 50 spins. 

III. SUMMARY 

In summary, we presented ESR spectroscopy using a single NV center in diamond. First, 
we demonstrated the identification of microscopic spin baths surrounding a single NV center 
and the investigation of static and dynamic properties of the bath spins using Rabi, FID, SE 
measurements as well as NV-based ESR spectroscopy. We also performed the investigation 
with several other single NV centers in the diamond sample and found that the properties 
of the bath spins are unique to the NV centers. Finally, by analyzing the intensity of the 
NV-based ESR signal using the computer simulation, we estimated the detected spins in the 
DEER measurement to be < 50 spins. 
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